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Abstract 

We introduce the concept of stationary metastable states (SMS's), 
and give a prescription to study it using a restricted partition function 
formalism. This requires introducing a continuous entropy function S{E) 
even for a finite system, a standard practice in the literature though never 
clearly stated, so that it can be differentiated. The formalism ensures that 
SMS free energy exists all the way to T = 0, and remains stable. We intro- 
duce the concept of the reality condition, according to which the entropy 
S{T) of a set of coupled degrees of freedom must be non-negative. The 
entropy crisis, which does not affect stability, is identified as the viola- 
tion of the reality condition. We identify and validate rigorously, using 
general thermodynamic arguments, the following general thermodynamic 
mechanism behind the entropy crisis in SMS. The free energy -Fdis(T) of 
any SMS must be equal to the T — crystal free energy Eo at two dif- 
ferent temperatures T — 0, and T — Tcq > 0. Thus, the stability requires 
FdisiT) to possess a maximum at an intermediate but a strictly positive 
temperature Tk, where the energy is E = Ek- The SMS branch below 
Tk gives the entropy crisis and must be replaced by hand by an ideal 
glass free energy of constant energy Ek, and vanishing entropy. Hence, 
Tk > represents the Kauzmann temperature. The ideal glass energy Ek 
is higher than the crystal energy Eo at absolute zero, which is in agree- 
ment with the experimenatal fact that the extrapolated energy of a real 
glass at T = is higher than its T = crystal energy. We confirm the 
general predictions by two exact calculations, one of which is not mean- 
field. The calculations clearly show that the notion of SMS is not only 
not vaccuous, but also not a consequence of a mean-field analysis. They 
also show that certain folklore cannot be substantiated. 
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1 Introduction 



Glass transition is a ubiquitous plienomenon^ believed to occur in a su- 
percooled liquid (SCL), which is one of the metastable states in Nature, and 
has been investigated for at least over eight decades since the earliest works of 
Nernstl^j and Simon0]. Despite this, a complete understanding of the transi- 
tion itself and related issues is still far from complete, although major progress 
has been made recently [S] [3 Theoretical and experimental investiga- 
tions invariably require applying (time-independent) thermodynamics to SCL 
to extract quantities like the entropy. This presumes, as we do here, that there 
exists a stationary limit of the metastable SCL state (SMS) under infinitely slow 
cooling in which the crystal (CR) is forbidden from nucleating 9 . We further 
assume that the SCL free energy can be defined, at least mathematically (see 
below for details), all the way down to absolute zero, which may not always 
be possible. It was discovered recently that under some conditions, SCL free 
energy can terminate in a spinodal at a non-zero temperature as we lower the 
temperature 

Experimentally, the SCL configurational entropy exhibits a rapid drop near 
a temperature ~ two-thirds of the melting temperature TmI^ El- K is found 
that the smooth low-temperature extrapolation of the measured excess entropy 
Scx{T) = S'scl(T) — Scr{T) becomes negative P below a non-zero temper- 
ature. Since it is hard to imagine CR, being more ordered, having a larger 
entropy than SCL, Kauzmann suggested that something like a glass transi- 
tion must intervene to avoid this entropy crisis, known commonly as the Kauz- 
mann paradox (S'cx(r) < 0)^, at some positive temperature. There are several 
computationaljSj and theoretical[7| 1,8, IQ, 12j| results clearly demonstrating the 
existence of some kind of entropy crisis. 

It should be stressed that there is no thermodynamic requirement for 5'cx(T) 
to be non-negative. There are physical systems like He^ in which S'ox(T') can 
become negative at low temperatures. A recent exact calculation on a classical 
system|3 clearly demonstrates a negative Scx{T) at low temperatures. If there 
is any hope of finding a thermodynamic basis for the glass transition, we must 
look for a condition for the glass transition that is dictated by thermodynam- 
ics. Thus, in the following, we replace Sc^{T) by the entropy S{T), such as 
the configurational entropy, which represents the natural log of the number of 
microstates W{E), where E is the average energy at that temperature. Conse- 
quently, a state with negative entropy is impossible to observe in Nature. The 
violation of the reality condition S{T) > signifies a genuine or absolute en- 
tropy crisis in Nature. We will interpret the entropy crisis in this work to signify 
the reality condition violation S{T) < 0, and denote the temperature by Tk, 
the Kauzmann temperature, when the violation begins as the temperature is 
reduced. 

It is widely recognized that Tk is a theoretical point and not accessible by 
experiment. However, its accessibility in itself is not important if its usage helps 
us understand or explain glassy behavior. It is well known that absolute zero is 
inaccessible; yet the study of a statistical model at T = is a time-honored first 
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step to study the physics of the model at higher temperatures. The concept 
of a Kauzmann temperature enables us to explain many glassy behavior. For 
example, the existence of the entropy crisis at a non-zero temperature is con- 
ceptually necessary for the observed Vogel-Fulcher behavior in fragile systems. 
Its importance cannot be denied just because it is experimentally inaccessible. 
A majority of the experts in the field continues to find the Kauzmann temper- 
ature to be an extremely useful concept. The most important consequence of a 
positive Tk is that there is a lower limit to supercooling before an experimental 
glass transition must intervene, and this lower limit is not absolute zero. This 
is an extremely useful information for experimentalists who are interested in 
investigating glassy behavior. 

The following two experimental observations play a very important role in 
our understanding of the glass transition: 

(Exptl) The heat capacity of the glassy state can be substantially different 
from that of the corresponding crystal at the same temperature |14|. 

(Expt2) The energy of the glassy state (Ek) after extrapolation to absolute 
zero is higher than that of the CR energy -EqD HI- 

Based on these observations, we have recently proposed |71 |H| a mechanism 
behind the entropy crisis and the ideal glass transition in the SMS of a system 
of long polymers. The proposed mechanism is as follows. 

Entropy Crisis Mechanism Since Ey^ > Eq, the SMS free energy rises 
rapidly below the melting temperature Tm as the system is cooled and crosses 
the {T = 0)—CR free energy Eg at some positive temperature 0<Tcq < Tm- The 
SMS free energy again equals Eq at T = 0. Therefore, it must go through a 
maximum at 0<Tk < Teq, at which the entropy of the system vanishes. Since 
negative entropy below Tk is not possible for physical states, the SMS below Tk 
is replaced by an ideal glass with a constant free energy Ek, and zero entropy 
and heat capacity. 

Since the above two experimental observations are the main ingredients for 
the proposed mechanism, it appears likely that the mechanism is common to 
all glass forming systems with metastability, and not just restricted to long 
polymers. If true, this will suggest that the above mechanism is a generic root 
cause of the entropy crisis at a positive temperature in all systems in which 
glasses are formed by cooling their metastable states. 

However, the situation has become very confusing, as there have appeared 
several recent arguments [TKl [T7I IH^ against the existence of the entropy cri- 
sis. These arguments contradict many exact calculations [3 |H1 1101 1111 El and 
simulations [B| that clearly establish a positive-temperature entropy crisis. We 
should also mention classical real gases, like the van der Waals gas, or the classi- 
cal ideal gas that are known to give entropy crisis at positive temperatures [T5] . 
It is, therefore, extremely important to clarify the issue, which we do by prov- 
ing the general validity of the above mechanism behind all glass formation in 
metastable states. Two important aspects in the mechanism require general 
justification 

(Ml) The SMS free energy, obtained by mathematical continuation of the 
disordered equilibrium state below the melting temperature, and the CR free 
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energy are the same at absolute zero. 

(M2) The temperature Tk when the entropy crisis occurs first is non-zero. 

Together, they support the existence of an entropy crisis at a positive tem- 
perature. The importance of a mathematical proof is that it settles the issue 
once for all. However, another very important consequence of the proof is, as we 
will see here, that if a phase in some calculation gives rise to an entropy crisis, 
it means that there must be another (ordered) phase without an entropy crisis. 
The proof uses general but rigorous thermodynamic and statistical mechanical 
arguments, valid for classical or quantum systems. We verify our conclusions 
by two exact calculations, one of which is not a mean-field calculation. 

As the issue has remained unresolved for so long most certainly implies that 
its resolution is not trivial. Thus, it should come as no surprise to the reader that 
our arguments are somewhat involved, can be divided into several clear parts, 
and require patience to go through. However, we believe that they are certainly 
not beyond the reach of a majority of the workers studying glasses. The proof 
of the validity of the above mechanism follows from the various theorems that 
we prove here. To guide the reader, we summarize the strategy for the eventual 
demonstration. 

1. We first prove that for > Eq, where i?D is the lowest energy for which 
SMS exists [S{E) > 0], the SMS free energy must become equal to the (T = 0)- 
CR free energy Eq at a positive temperature Teq. Neither the SMS entropy nor 
its inverse slope at Eu need vanish for this to hold. 

2. We then prove that if the energy of the SMS state is higher than Eq, then 
its temperature must be positive. From this it follows that SMS cannot reach 
T = with an energy E — En > Eq . 

3. We then prove that SMS free energy at T = again equals Eq at T = 0. 
Thus, the SMS free energy must have a maximum somewhere in the range 
< T < Teq. The entropy at the maximum is zero, identifying the maximum 
with the Kauzmann point and its temperature with Tk. This also identifies the 
lowest energy Ey) with T^Ki the energy at Tk. This energy is greater than Eq 
because of the non-negative SMS heat capacity. 

The existence of the maximum in the SMS free energy thus finally proves 
the existence of a Kauzmann point. 

1.1 Fundamental Postulate 

Our general proof assumes the existence of SMS's, so that thermodynamics can 
be applied. The need for the assumption is easy to understand. At present, 
our understanding of whether equilibrium states can be demonstrated to exist 
mathematically even in simple models is too limited. We should recall that the 
existence of equilibrium states is taken for granted as a postulate in statistical 
mechanics and thermodynamics, where it is well known that it is extremely hard 
to prove their existence. We quote Huang[2D]: "Statistical mechanics, however, 
does not describe how a system approaches equilibrium, nor does it determine 
whether a system can ever be found to be in equilibrium. It merely states what 
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the equilibrium situation is for a given system." RuellejU notes that equihb- 
rium states are defined operationally by assuming that the state of an isolated 
system tends to an equilibrium state as time tends to +00. Whether a real 
system actually approaches this state cannot be answered. We make a similar 
assumption about the existence of SMS within the restricted ensemble that will 
be introduced below. Any analysis of thermal data on glasses and supercooled 
liquids requires using time-independent thermodynamics, as discussed above. 
Thus, hypothesizing the existence of SMS is perfectly justified by the practice 
in the field. The existence of supercooled liquids, and glasses in systems with 
short- or long-range interactions or structural glasses, and their stability over a 
long periods of time is undeniable. These stable states are a manifestation of 
SMS in these systems. 

Without such an assumption, wc cannot justify using conventional (time- 
independent) thermodynamics to analyze SCL data. The two exact calculations 
in this work show that the hypothesis is not vacuous. The behavior in real 
systems is, of course, oblivious to the state of our knowledge and is not controlled 
by the fact that we can only demonstrate SMS in a few exact calculations. 

1.2 Reality Condition and Entropy Crisis 

As said above, the entropy S(T), such as the configurational entropy, that we 
consider in this work represents the natural log of the number of microstates 
W{E) where E is the average energy at that temperature. For the microstates 
to exist in Nature, it is evident that W{E) must satisfy the reality condi- 
tion W{E) > 1 [S(T) > 0]. However, a state with negative entropy can emerge 
in theoretical calculations or extrapolations. If it happens that the calculations 
or extrapolations result in a negative S(T) below some positive temperature 
Tk, this should be interpreted as the absence of real microstates available to 
the system at those temperatures and the system cannot be found in those mi- 
crostates in Nature. It is this violation of the reality condition S{T) > that 
signifies the genuine entropy crisis in Nature. 

It should be noted that our criterion for the entropy crisis is much more 
stringent than the original Kauzmann requirement ScxiT) < 0, since it is pos- 
sible to have S{T) > even though Scx{T) < 0. Thus, our Tk is lower than the 
corresponding temperature where ScxiT) — 0. 

There are two independent aspects of thermodynamics and statistical me- 
chanics. The first one is the requirement of stability according to which thermo- 
dynamic quantities like the heat capacity, the compressibility, etc. must never 
be negative. The other aspect, independent of the stability criteria, is the reality 
condition that ensures that such states occur in Nature[221- The reality issue is 
central in our approach and is discussed further below. We will see below that 
the mathematical extension of the free energy of the disordered liquid phase 
below the melting temperature Tm will always satisfy the stability criteria ev- 
erywhere (T > 0), but the reality condition is satisfied only above some positive 
temperature Tk- 
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1.3 Common Folklore 

Majority of the calculations displaying stationary metastable states are in spin 
models, and are carried out at the level of the random-mixing approximation 
(RMA); the latter is valid in the limit of infinite coordination number q and 
vanishing interaction strength J, keeping qj is fixed and finite. This approx- 
imation is equivalent to solving the models in an infinite- dimensional space. 
This has given rise to the common folklore that SMS's occur only in an infinite- 
dimensional space whose coordination number is also infinite. This is incorrect 
as we will clearly demonstrate here by an explicit calculation. The calculation 
is a non-mean-field calculation carried out exactly in a 1-dimensional model 
and captures SMS. The calculation is presented to overcome this folklore. The 
model also shows the entropy crisis in SMS. Another alternative way to visual- 
ize the RMA is to think of long-range interactions in a finite-dimensional space. 
Thus, another folklore is that SMS does not exist for short range models. Even 
frustration is considered in the folklore to be necessary for the glassy behav- 
ior. The second example in the work is that of an Ising model with short-range 
interactions and no frustration. This example establishes not only the exis- 
tence of SMS in short range models but also that frustration is not necessary 
for the existence. The model is solved exactly on a Husimi cactus; hence, its 
thermodynamics is proper. 

It is true that the cactus can only be embedded in an infinite dimensional 
space, but locally it resembles a square lattice which has a finite coordination 
number. The important point to note is that the interactions are short-ranged. 
The dimensionality is relevant only if we are interested in critical exponents, 
which we are not in this work. The finite coordination number of the cac- 
tus makes our calculation very different from RMA. Our general proofs are 
certainly not based on RMA ideas. From all the experience we have accu- 
mulated, the cactus provides a much better description of the square lattice 
model than the conventional mean-field (RMA) , as we have shown elsewhere. [23 



2 Formulation 

2.1 Canonical Partition Function 

The stationary nature of the SMS allows us to investigate it using the partition 
function (PF) formalism. We consider a system composed of N particles con- 
fined in a given volume V and at a given temperature T. The thermodynamic 
limit is obtained by taking the simultaneous limits N oo, and V oo, such 
that the volume per particle v = V/N is kept fixed and finite. The internal 
degrees of freedom of the system contain the configurational (i.e., positional) 
degrees of freedom and all other degrees of freedom like the translational de- 
grees of freedom coupled to them[23. The canonical PF determined by these 
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degrees of freedom is given by 

Zn(T) = Tr eM-PE), 



(1) 



where the trace operation Tr is over the coupled degrees of freedom and /? = 1/T, 
T being the system temperature in the units of the Boltzmann constant fee. In 
the following, we will usually suppress the index N on the PF, unless necessary. 
The energy of the system E is determined by the coupled degrees of freedom. 
The PF in is irreducible in that it cannot be written as a product of sev- 
eral non-trivial PF's corresponding to independent sets of degrees of freedom; 
see 123. Some coupling between different degrees of freedom, no matter how 
weak, is required to achieve equilibrium so that they all share the same common 
temperature. The weaker the coupling, the longer the time required to come 
to equilibrium. The temperatures in different PF factors in the product need 
not be the same as there is no coupling between independent degrees of free- 
dom. Therefore, such a situation does not have to be considered here. We only 
consider the case of an irreducible PF in this work. 

The microstates over which the operation Tr is carried out in |^ are deter- 
mined by the coupled degrees of freedom, while the macrostates are determined 
by the thermodynamic variables N, V and the temperature T. Since the degrees 
of freedom are not part of the macroscopic description, but the average energy 
defined below is, it is convenient to replace the trace operation in J^l by a trace 
only over the energy of each microstates. We classify different realizations of 
the degrees of freedom, i.e. the microstates, according to their energy so 
that W{E) is the number of microstates of energy E. We can now transform 
the above PF in ^ into 

Z(T) = Tr W{E) exp(-/3£:), (2) 

where Tr now stands for the trace operation over the energy (eigenvalues) E > 
Eq up to its maximum. We also introduce the following extensive quantities 

n{T) = \nZ, F{T) = -T\nZ; (3) 

here, F(T) is the conventional Helmholtz free energy. 

Since the sign of the entropy is an important issue in the study of glasses, 
it is important that the entropy be introduced using the number of microstates 
W{E), so that W{E) > 1. This requires some kind of discretization of the 
degrees of freedom PSI to count the microstates. In the following, we will 
assume that such a discretization has been carried out so that we always have 
W{E) > 1 for states that occur in Nature. This ensures that the corresponding 
microstate entropy S{E) = \xiW{E) due to the coupled degrees of freedom is 
non-negative. Despite this, we will see below that negative entropies can be 
obtained from the free energy associated with the SMS at low temperatures, 
even though the free energy itself is stable ■ This will force us to conclude 
that an ideal glass transition must intervene to avoid the resulting entropy crisis. 

The PF in (0) is irreducible. It can happen in some cases that the PF is a 
product of several non-trivial PF's corresponding to independent sets of degrees 
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of freedom, each set containing only coupled degrees of freedom; see |23- The 
application of the reality condition to each set requires that the entropy from 
each set be non-negative for the system to occur in Nature. The entropy crisis 
occurs when any of the entropies from the independent sets (each containing 
coupled degrees of freedom) becomes negative even if the total entropy due to 
all sets is non-negative. 

2.2 Thermodynamic Limit 

The thermodynamic limit is obtained by taking N oo, keeping v = V/N 
fixed. To make the discussion clear, we will exhibit the subscript N in various 
quantities in this subsection. The limit is taken by considering the sequence 
formed by 

ujn{T) = {l/N)nN{T), where n^iT) = IhZnIT), 

for different values of iV as iV ^ cx). The volume may be changed according 
to V = vN. For proper thermodynamics, the limit of the sequence must exist, 
which we denote by to. 

We express the fact that V/N is fixed by stating that F is a homogeneous 
function of order 1 in N. Similarly, the existence of the limit w is expressed by 
stating that ilN{T) or Fjv(T) = —TIuZm are homogeneous functions of order 
1 in N. One can also say that the temperature T is a homogeneous function of 
order in N. The average energy Em{T) and the entropy Sn{T) are homoge- 
neous functions of order 1 in N. For any quantity Qn, which is a homogeneous 
function of order 1 in TV, the ratio Q n /N is a homogeneous function of order 
in N and must possess a limit q as ^ oo. We express this fact in the following 
as 

Qn/N ~ q, (4) 

whose significance is as follows: for finite N, the meaning of is that the right 
side may differ from q, but the difference vanishes as iV ^ cxd. 

Remark 1: In the following, whenever we compare different extensive quan- 
tities Qi or different PF's Zi, it is implicit that we are comparing the ratios Qi/N 
or hiZi/N, respectively. 

2.3 Continuous Energy and Entropy Functions for Finite 

The average energy and entropy in the canonical ensemble are continuous func- 
tions of T and are given by 

E{T) EE ~{dnid(3), (5a) 
S{T) EE ~{dF/dT), (5b) 

respectively; we have suppressed the subscript N for simplicity. They should 
not be confused with the microstate energy E and entropy S{E). To see this 
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clearly, we rewrite (j5b|l as S = In Z + E/T, and exponentiate it. We can use the 
explicit T-dependence of E(T) to express S{T) as an explicit function S{E) of 
E = E{T). Introducing W(E) = exp[5(!E)], we have 

W(E) ^ Tr W{E)e-l^^^~^^ = W{E,r,)e-l^^^"^-^'> + Tr' W{E)e->^'^^~~^\ (6) 

where the prime on the trace indicates that some microstate energy E — E-^^ is 
not allowed in the trace. For a finite system {N < oo), the allowed microstate 
energy E, and the microstate entropy S{E) = \nW{E) are discrete, while the 
average energy E and the average entropy S are continuous functions. Let us 
consider the case when the average energy is exactly equal to the microstate 
energy E^. From lO, we observe that 

W(E) > W{E^). (7) 

The left side in Q represents the value of the continuous function W{E = -Em), 
while the right side is the discrete quantity W{E,n) for the finite system. The 
difference between them is due to the last term in which is non- negative, 
which vanishes as N —f oo. For finite TV, the microstate entropy S{E) represent 
isolated points, which get closer and closer to the continuous function S{E) as 
N increases, and cover it entirely as N oo. Similarly, E represents a contin- 
uous variable, and contains in its range isolated values given by the microstate 
energies; the latter cover the entire range of £' as — ^ oo . 

The continuous function S{E) is amenable to differentiation, which is not 
possible with the discrete set of points S{E), and contains all the useful physical 
information. Because of this, there is no harm in expressing S{E) simply by 
S{E)] the latter now represents a continuous function. Similarly, there is no 
harm in expressing the continuous variable E by E. This is a common practice 
in the literature. Whenever we need to make a distinction between the discrete 
values and the continuous functions, we will speak of the values related to mi- 
crostates in the former case, and of average values in the latter case. Similarly, 
we use S{T) to express the entropy S{T) = S{E) as a function of T, where 
T is the temperature at which the average energy from (|5a|l is E. The use of 
the continuous functions enables us to consider the entropy for any value of the 
energy, whether it represents the energy of a microstate or not. The continuous 
entropy function satisfies the following relation: 

(dS/dE) = {dS/dT)/{dE/dT) = 1 /T. (8) 

Let us first illustrate the above points by a simple example. Consider two fer- 
romagnetically interacting (interaction strength — J) Ising spins in the absence 
of any external magnetic field. There are four microstates, two with E = — J 
(the lowest possible energy), and two with E — J (the highest possible energy). 
Thus, there are only two microstate energies, each having the microstate entropy 
In 2. On the other hand, the average energy and entropy are continuous func- 
tions of T. To see this, we calculate the canonical PF, which \s Z = 4 cosh if. 
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where K = J /T. A simple calculation for the average energy and entropy yields 

E = - J tanh K, S = hi(4 cosh K) - K tanh K. 

We observe that the energy ranges continuously between — J and , and the 
entropy ranges continuously between In 2 and In 4 for T > 0. (The energy range 
between and J corresponds to negative temperatures that we do not consider 
here.) The average energy and entropy are (—J) and In 2, respectively, as T ^ 0, 
and and In 4, respectively, as T — > oo, as expected. The average energy, and 
entropy fill continuously the microstate energy and entropy gap, respectively. 
The important point is that the lowest energy corresponds to the absolute zero 
of the temperature. 

We now prove that the entropy slope relation in ^ is valid for any N , finite 
or not. We consider the PF in Q for finite N. We start with the continuous 
entropy and energy functions related by TS ^ E — F. Differentiating on both 
sides with respect to T, we find 

S + T{dS/dT) = {dE/dT) + lnZ + T{dZ/dT)/Z. 

Recognizing that 

(dZ/dT) = (l/T2)Tr EW{E) eii-p{-f3E), 

we find that the above equation reduces to 

T{dS/dT) = (dE/dT), 

which proves © for finite N. The limit ^ oo is considered by dividing both 
sides of TS ~ E — F hy N, and carrying out the same steps. 

We now argue that if we replace W{E) by W{E) in (O to define a new PF 

Z{T) = Tr W(E)e~'^'^, (9) 

in which the trace is still over the microstate energies, then using S = In Z+E/T, 
and W{E) = exp(S'), we find that 

Z{T) = Z{T)Tr 1 ~ Z(T), (10) 

since Tr 1 is the number of distinct microstate energies and satisfies ln(Tr l)/iV ~ 
0p5|. Consequently, lnZ(T)/A^ and \nZ{T)/N are the same in the sense of the 
above Remark 1. In the following, we will consider both versions of the PF for 
macroscopically large but finite systems, as they are identical in all thermody- 
namic consequences. Indeed, it is a common practice in the literature not to 
make any distinction at all. Because of this, we will denote both of them by 
Z(T); we will indicate the difference whenever needed. 

Remark 2: It should be noted that the association of the slope in JSJ 
with the temperature T requires fixing the temperature of the system. The 
temperature is fixed from outside, such as by using a heat bath, and must be 
independent of the size N of the system, even though both S and E depend on 
N in accordance with 10}. In this sense, the significance of T = l/{dS/dE) is 
that for any N, there exists an E so that ((HJ is always satisfied. 
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2.4 Configurational PF 

In classical statistical mechanics (CSM), the positional degrees of freedom are 
independent of the translational (i.e., momentum) degrees of freedom when the 
collisions are neglected. In this case, the PF can be written as a product of 
two PF's. One of them is ^keCTke) determined by the momentum degrees of 
freedom; here Tke is the temperature associated with these degrees of freedom 
and the energy associated with these degrees of freedom is the kinetic energy of 
the system. (It should be noted that the momentum degrees of freedom are also 
independent of each other so each momentum degree of freedom can have its 
own temperature. We will not worry about this complication.) The other factor 
is the PF Z{T), known conventionally as the configurational PF at a tempera- 
ture T. It is determined only by the positional degrees of freedom for which the 
energy is the potential energy in the system. The collisions between particles 
are neglected so the two sets of the degrees of freedom have no mechanism to 
come to equilibrium. Hence, there is no way for the two temperatures to be 
the same. Usually, it is assumed that there exists a weak coupling between the 
two degrees of freedom, which is sufficient to bring about eventual equilibrium 
between them. Thus, to a good approximation, the above factorization is taken 
to be valid in real systems. The important observation is that the value of Zke 
is independent of the value of the configurational PF Z(T) at this level of the 
approximation. Consequently, the total entropy Stot{T) due to both energies is 
additive: S'tot(r) = S{T) + S'KE(r), where S'ke(T) is the entropy due to the 
kinetic energy and is independent of the configurations of the system. Further- 
more, Ske{T) is the same for all kinds of systems. Thus, Zke{T) is not of any 
interest when studying any particular system. Because of this, there is no harm 
in restricting our attention to the studying the configurational PF Z{T). In this 
case, the PF in (0) represents the configurational PF Z{T) as defined conven- 
tionally in CSM so that S{T) will represent the conventional configurational 
entropy in CSM. The entropy crisis occurs when S{T) becomes negative, even 
if Stot(T) is non- negative. 

2.5 Quantum PF 

In quantum statistical mechanics (QSM), the kinetic energy is an operator and 
cannot be separated out from the total PF|25. The role of the configurational 
PF is now played by the total PF. The energy E va ^ now represents the 
eigenvalues of the total (potential+ kinetic) energy, and S{T) derived from ((2Jl 
now represents the total entropy, which cannot be broken into additive terms as 
was the case in CSM above. Thus, one cannot define the classical configurational 
entropy in QSM. We can think of S{T) derived from as the quantum analog 
of this classical concept. 

The irreducible PF Z{T) in ||2Jl is the general form of the PF valid in both 
CSM and QSM, with S{T) = S{E) equal to the conventional configurational 
entropy in CSM, and the total entropy in QSM. From now onward, we will 
no longer explicitly distinguish between the classical and quantum PF's. Our 
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discussion is going to be valid for both cases. 

2.6 Conditions for Equilibrium and Negative Entropy 

For the lowest allowed energy Eq, we must surely have W{Eq) ^ 0. Assuming 
TS{T) — > as T ^ 0, we recognize that Eq represents not only the Helmholtz 
free energy but also the energy of the perfect CR at T = 0. (We assume that CR 
has the lowest free energy at low temperatures.) Since W{E) is non-negative, Z 
is a sum of positive terms. Because of this, the probability of every microstate is 
strictly non-negative. As a consequence, the following two principles are always 
satisfied. 

(1) Maximization Principle. The PF Z must be maximized in the ther- 
modynamic limit N — > GO, keeping V/N fixed. The maximum value of Z(T) 
corresponds to picking out the maximum term e^~^^ in |(2Jl. This maximum 
term corresponds to E = E, defined above. 

(2) Stability Principle. The heat capacity, which is given by thy fluctuations 
in the energy is non-negative. This remains true even in the thermodynamic 
limit. 

It should be stressed that the non-negativity of the heat capacity and the 
maximization principle only require the positivity of W{E). Thus, both princi- 
ples remain valid even if the entropy becomes negative j22j. Stability and reality 
are two independent aspects of our formalism. This observation is going to be 
useful when we discuss the metastable states below. 

3 Stationary Metastable States and Restricted 
Ensemble 

3.1 Infinite System 

Conventional statistical mechanics describes equilibrium states, which satisfy 
the above two principles of reality and global maximization. For this, it is 
necessary that we have N oo. The existence of a melting transition, which 
also requires N ^ oo, at Tm means that the disordered equilibrium liquid (EL) 
phase above Tm and the ordered CR below Tm correspond to different values of 
the order parameter p, which is traditionally used to distinguish various phases 
of the system, with p — representing the disordered phase and p ^ the 
ordered phase CR. (One of our examples below will show explicitly how the 
microstates can be divided into the two disjoint classes.) We assume here for 
simplicity that there is only one kind of ordered phase. The extension to the 
case of many disjoint ordered states with different non-zero values of p is easy 
to incorporate in the approach. This distinction in the order parameter is easily 
made in the case iV ^ oo, by denoting the free energy per particle {N oo) 
above Tm by /dis(T), and below Tm by /ord(T), from which we can calculate 
the entropies, and energies per particle 

Sa{T) = ^{dfc/dT), CciT) = -{d(3U/dl3), a = dis, ord, (11) 
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respectively, corresponding to the two states. Due to the limit N ^ oo, the 
per particle microstate entropies and continuous entropy functions are equal at 
the same e. Thus, we can think of the continuous entropy and energy functions 
Sct(-E), and ea{T) for the two states. 

The global maximization, which operates when N oo, is required to argue 
that /ord(r) is the equilibrium state free energy (for CR) below Tm, and fdis{T) 
the equilibrium state free energy (for EL) above Tm : 



Thus, incorporating the global maximization principle will not allow us to de- 
scribe SCL. The singularity in the equilibrium free energy per particle f(T) at 
Tm in (|12all2b|l . which forces it to switch from fdis{T) to foTd{T) at Tm, is a 
hallmark of the a phase transition. This singularity in /(T) does not necessarily 
imply a singularity in either of its two pieces /dis(T') and ford{T). Both of them 
can exist on either side of Tm. This possible extension is not a consequence of a 
mean-field approximation, as our first example will demonstrate. A consequence 
of this is that the functions Sa{T), and ea(T) also exist on either side of Tm. In 
the following, we are only interested in the case in which Sdis(T), and edis(T) 
exist all the way down to T = 0. 

A prescription to describe metastability using the PF formalism can now be 
formulated. 

Metastability Prescription We abandon the global maximization princi- 
ple, and use fdis{T) to give the free energy of the metastable disordered phase 
(supercooled liquid) below Tm and /oid(r) to give the metastable (superheated 
crystal) state free energy above Tm. Similarly, Sdis(r), edis(T) and SordiT), 
edis(T') give the entropy and energy per particle for the supercooled liquid and 
superheated crystal, respectively. 

However, in this work, we are only interested in the supercooled liquid. 

3.2 Finite System 

The study of an infinite system allows us to identify ordered (p ^ 0) and dis- 
ordered (p = 0) microstates and gives their entropy and energy per particle 
Sc((T), and ea{T), respectively. The identification is useful to classify each mi- 
crostate as ordered (p ^ 0) or disordered (p = 0). This does not imply that 
it is feasible to count these microstates for an infinite system in which there 
are infinite microstates. For a complete analysis, we need to be able to count 
in principle the microstates and classify them. For this, we need to consider 
the case of finite but large N. In this case, the continuous entropy and energy 
functions for the entire system are SaiT) ~ A^Sq(T), and Ea{T) ~ NcaiT). We 
also consider the continuous functions Sa{E), and Wa{E) = exp[Sa{E)] which, 
as we show later, exist for all energies E > Eq. As said earlier, the continuous 
function Wa (E) for microstate energy E is close to the number of ordered and 
disordered microstates of microstate energy E; see {T)). 



/(T) = /o,d(T) T<Tm, 
/(T) = /dis(T) T>Tm. 



(12a) 
(12b) 
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The total number of microstates W{E) with microstate energy E can be 
written as a sum of the number of microstates Wdis{E) consistent with p = 0, 
and the number of microstates WoidiE) consistent with p 7^ 0, so that 



W{E) = Word 



{E) + Wdis{E). 



(13) 



(Presence of other configurations will not affect our argument.) While Woi-d{E) 
certainly exists for microstate energies starting from E = Eq, there is no 
guarantee that VFdis(-E) also exists near E — Eq. Most probably, microstate 
Wdis{E) > 1 does not continue all the way down to E = Eq. If it did, the 
energy of the disordered phase at absolute zero would be Eq, the same as that 
of CR. This would most certainly imply that they would coexist at T = 0, each 
having the same volume; recall that we are considering a fixed volume ensemble. 
While there is no thermodynamic argument against it, it does not seem to be 
the case normally. Usually, the most stable state at T = is that of a crystal. 
Moreover, it is an experimental factH^ that all glasses have much higher ener- 
gies or enthalpies compared to their crystalline forms at low temperatures; see 
(Expt2) above. Anticipating that the ideal glass, which is the stationary limit 
of all observed glasses, should have its energy higher than E'o, we conclude that 
the lowest possible energy i?K for the disordered state is strictly larger than Eq. 
In other words, the microstate number Wdis{E) has the following property: 



On the other hand, we are only interested in considering the case when Sdis(r) 
exists for ah T > 0. Thus, S'dis(r) exists for aU T > 0. We will show later 
that the temperature T = Tk — l/{dSdis/dE) at Ek > Eq corresponding to 
the disordered phase is strictly positive. Thus, Sdis{T) [or S'dis(i?)] must be 
negative below T — Tk [or E = Ek]- We will, in fact, see that Sdis{E) continues 
all the way down to Eq. 

3.3 Restricted Ensemble 

We continue with finite but large N. Using the continuous forms of Woid{E) and 
WdisiE), we introduce the following restricted ensemble approach to describe 
metastable states^j. We follow (O and we introduce two new PF's using by 
replacing continuous W{E) by continuous Word{E) and Wdis{E), respectively: 



WdisiE) =0 for E<Ek, 
Wdis{E) > 1 for E>Ek. 



(14a) 
(14b) 



Zc{T) EE TrWa{E) exp{-PE), a = dis, ord. 



(15) 



and the corresponding free energies 



T\i\Za{T), a = dis, ord. 



(16) 



As A^->(X3, F„(T)/7V^/„(T). 
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Above the melting temperature Tm, .^dis >^ord, and below Tm, ^dis<-^ord- 
The use of the maximization principle describes CR and EL collectively for 
thermodynamically large N: 

Zdis(T) > Zo,.d(T), T > Tm, (17a) 

Zord(T) > Zdi«(T), T < Tm, (17b) 

so that 

Ford{T) > Fdis(T), T > Tm, (17c) 

Tdis(T) > To,d(r), T<Tm. (17d) 

For finite A^, there cannot be any singularity at Tm, so that Tdis(T) or 
^ord(r) has no singularity of its own at Tm, i.e., they exist on both sides of Tm- 
We now describe SMS by the PF Zdis(T) below Tm - The continuation necessarily 
yields p = also for SCL. 

The entropy functions Sa{E) are shown in Fig. 1: OHAB represents S'ord(T'), 
and KCO AD represents 5dis(T')- The entropy as a function of E must be 
thought of as the entropy in the microcanonical ensemble 25 , which must be 
maximum for the equilibrium state. Since a SMS is not an equilibrium state in 
the unrestricted ensemble, its entropy at some E cannot exceed the entropy of 
the corresponding equilibrium state at the same E. A consequence of this en- 
tropy condition is that the free energy Tdis(T) of SMS cannot be lower than the 
free energy Toid (T) of CR at the same temperature T. This explains the form of 
the entropy and free energy in Fig. 1 . The two free energies are sown in the inset 
by OAB representing Tord(7"), and OKCO D representing Fdis(T). The form of 
the entropy and free energy is also supported by all known observations^] 
exact 13 IHl 112) and numerical"^ calculations, and from the arguments and the 
calculations presented below. We note that 

SordiE) < SdisiE), E > Em, (18a) 
Sdis{E) < Sord{E), E < Em, (18b) 

where Em is the energy at A where Sord{E) = Sdis{E); see Fig. 1 . The SMS 
corresponding to the stationary SCL is defined by the branch KCAH . Similarly, 
superheated CR is defined by the branch HAB. We note that the entropy 5dis 
of the metastable branch goes to zero at Tk>0. This behavior is supported 
by our rigorous proof and by two exact calculations in the paper. However, 
we will allow the possibility in our discussion below that the lowest energy for 
Sdis is Ed, and that Sdis{Eu) > 0. We also show the entropy function for a 
non-stationary metastable state by GF in Fig. 1, assuming that crystallization 
is forbidden. Wc do not consider the non-stationary state GF anymore. 

3.4 A Useful Theorem 

We now prove an extremely useful theorem. 
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Theorem 1 Since Ek > E^, the free energy Fdis{Teq) at O equals the free 
energy Ford{T = 0) = Eq at O, where T^q is the inverse of the slope of the 
line 00 touching the entropy function Sdis{E), which vanishes at E = Ek; see 
point O in the inset in Fig. 1. 

The proof is very simple. The slope l/Toq of 00 is given by 

l/Teq = ^dis(£;o')/(^^0' -^O), 

where Eqi is the energy at O . Thus, 

= EqI — T'cqS'dis(£'Q' ). 

Since the slope of 5dis(£'o') at Eq' is l/Toq, the right side represents the free 
energy Fdis(T'eq) of the SMS at O . The left side represents the free energy of 
CR at T = 0. This proves the theorem. 

It should be stressed that the proof does not use the vanishing of S'dis(£^K)- 
Thus, the equality -Fdis(7eq) = Eq is also valid if S'dis(£^D) > 0. The proof also 
does not depend on the entropy slope at E-o or i?K- The proof only requires 
that this slope be larger than or equal to 1/Toq. 

We now consider the behavior of Zdis(T'). Since the continuous entropy 
Sdm{E) is real, we have Wdis{E) > 0. Thus, we observe that Zdis{T) is a sum 
of positive terms. Because of this, Zdis{T) also satisfies the maximization and 
stability principles, just as Z{T) does. The only difference is that they are only 
valid in the restricted ensemble of the disordered microstates corresponding to 
P - 0. 

The reality condition cannot be violated, even for a SMS if it is observable. 
Therefore, its violation (S'dis(r) < for T < Tk > 0) implies that SMS cannot 
exist in Nature when the violation occurs, in which case the SMS associated 
with the SCL must be replaced by a new state, commonly known as the ideal 
glass state below Tk, whose energy at T < Tk > is Ek > i?oD I2- The 
transition between the two states is called the ideal glass transition. 

4 Finite Entropy Slope of the Disordered Phase 

Theorem 2 The slope of the entropy function Sdis at E/Eq > 1 must be finite. 
We again consider finite and infinite N separately. 

4.1 Finite 

We consider the average entropy functions Sa{E) for fixed V,N; see Fig. 1, 
where they are shown schematically as functions of E. According to (jHJ 

{dS/dE) = l/T, {dS^/dE) = l/T^, (19) 
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where 1/T, 1/Tq, represent the inverse temperatures corresponding to S, S'q,, 
respectively. By introducing the entropy and energy densities 

SN = S/N, s^,N = S^/N, CN = E/N, (20) 

we can rewrite (|19|l that will be useful in taking the thermodynamic limit below. 
For finite N, sn, Sa,Ar, and arc homogeneous functions of order 0; see I0J. 

The behavior of the slope in (|19|l for the disordered phase is critical in 
understanding what happens at K, the Kauzmann point (where S'dis = 0). There 
are two distinct possibilities. The slope at K is either finite, as shown explicitly 
in Fig. 1 and consistent with our rigorous analysis, or infinite. The former 
corresponds to a positive Tk, while the latter corresponds to an absence of a 
Kauzmann paradox (Tk = 0). Almost all explicit calculations show the former 
as the usual behavior, two of which are presented below. 

To proceed further, we consider (|13|l for any microstate energy. From what 
has been said earlier about the continuous function S{E), we extend H13(> to any 
allowed continuous energy so that it can be differentiated. Taking the derivative 
with respect to E on both sides, we find that 

W{dS/dE) = Word{dSo,d/dE) + Xd;,{dSdis/dE). 

We introduce the ratios = Wa/W (0 < < 1) and rewrite the above 
equation as 

1/T = Xoid I Tord + ^A\s /Tdis • (21) 

For T > Tm, ^ord ~ exp(S'ord - S'dis) << 1, while X^is ~ 1. For T < Tm, 
^dis ~ exp( S'dis — Sord) << 1, while Xord ~ 1. Let us apply this relation to 
the crystal at a temperature Tm > T > at which the average CR energy is 
E > Eq. It is clear that for any finite N, no matter how large, Xdis > 0. At the 
same time, the temperature T must be close to Tord of the CR, which is strictly 
positive, since E > Eq. This is possible only if 

Tdis>0 for E/Eo>l. (22) 

Thus, as long a,s E > Eq, the temperature of the disordered phase at E must 
be positive. This proves the theorem for finite N. 

4.2 Infinite 

Let us now take the limit iV — > oo. We only consider the interesting range 
2m > T > 0. The values of E and S keep changing with N for a given T. 
Thus, it is convenient to consider the sequences {ejv}, and {sq.at} for different 
N, so that {dsa.N Ide^) — 1/Tdis > 0. We now consider for larger and 
larger iV, such that Tdis > is kept fixed. From what was said earlier, it is 
clear that Xdis ^ from above, so that T — > To^d from below. But at every 
step of the hmit ^ oo, remains vahd. Thus, Tdis > for E/Eq > 1 
even when N oo. 
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The above proof of neither requires nor shows that the entropy Sdis{E) = 
for E > Eq. It also does not depend on the sign of Sdis{E) as long as it re- 
mains differentiable. It only requires Wa to be non-negative for Tq, to be real. 
Thus, it does not directly prove the existence of a Kauzmann point below which 
the entropy crisis would occur (Wdis < !)■ For this, we need to show that 
the entropy of the disordered phase becomes zero at some Ek > Eq. Then the 
above theorem proves that this occurs at a positive temperature. Furthermore, 
it shows that the disordered phase can reach the absolute zero only \i E = Eq. 
The proof that an entropy crisis does occur in general at a positive temperature 
is given in the following section. 

5 Free Energy at T = 

We now prove the following theorem. 

Theorem 3 The free energy Fa of all stable phases, mathematically continued 
or not, are equal in that F^/Eq —^ 1, provided Tsq. — > as T — > 0. Their 
entropies s^, however, may be different. 

The proof will require considering finite N and then taking the thermody- 
namic limit later just as above. 

5.1 Finite 

Consider finite N. From (|14a|l . we conclude that microscopic entropy Sdis{E) 
does not exist (as a bounded quantity) for E < Ek- On the other hand, it must 
be close to the continuous function Sdis{E) exists only for E > Ek, and vanishes 
at Ek, as shown in Fig. 1. This is consistent with (jl4b|) . From Theorem 2, we 
conclude that the temperature Tk corresponding to Sdis{E) at Ek is positive. 
Moreover, since N is finite, Sdis{E) is not singular at Ek. (Recall again that the 
singularity can only appear in the limit TV — > oo. Our examples will show that 
the entropy per site Sdis{E) remains non-singular even in the thermodynamic 
limit.) Then, it can be continued mathematically as a real function below Ek 
all the way down to Eq. But the continuation of S'dis(£') for E < Ek, which 
must necessarily lead to negative S'dis(£') there, is certainly not close to the 
microscopic entropy Sdis{E) {— —oo). Thus, the mathematical continuation 
of Sdis{E) below Ek will most certainly not represent the physics correctly. 
This is not surprising in view of the fact that it violates the reality condition. 
Despite this, the continuous entropy is still useful in the investigation, and 
we will allow it to become negative below Ek, as it allows us to draw a very 
important conclusion, as we demonstrate below. Corresponding, we will allow 
Wdis{E) < 1 under continuation. 

Due to the mathematical continuation of Wdis{E), both PF's Za{T), a —dis, 
ord, contain all energies from Eq upwards, so that both can be investigated in 
a unified fashion. As long as Wa{E) > 0, Za{T) is a sum of positive terms, so 
that ^aiT) = lnZa{T) satisfies stability principle. For example, the free energy 
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Fdis(T) is a concave function (a function which always lies above the line joining 
any two points on it) of T, as shown in the inset in Fig. 1 by OKCD. We take 
out the term corresponding to E = Eq from the trace operation, and express 

ZciT) = W^{Eo)e-^^«[l + Z'^{T)], (23) 

where we have introduced a new quantity 

Z'JT) = Tr'[W^iE)/Wa.{Eo)]e-^^''-''«\ (24) 

in which the trace operation Tr is restricted to all £^ > ii^o- It is assumed that 
WaiEo) > 0. Since E-Eq > 0, we note that e-^^^"^") as T 0. For finite 
N, Wa{E) is a bounded quantity, and so is the ratio Wa{E)/Wa{Eo). Hence, 
the product under the trace in H24|) vanishes, and so does Z^{T) as T ^ 0. We 
finally have 

Z„(r) ^ Wa{Eo)e~^^° as T ^ 0. 

We thus find that F„(r) Eq - TSa{Eo) as T 0. From the boundedness of 
Wa{Eo) > 0, we also conclude that TSa ^ as T — > 0; hence, 

F^/Eo ^1 as T (25) 

for CR and SMS. 
5.2 Infinite 

We recognize that both Fa and Eq are homogeneous functions of order 1 in 
N. Thus, their ratio Fa/Eo is a homogeneous functions of order in N. Conse- 
quently, we can take the thermodynamic limit A'^ — > cxd by exploiting Q without 
altering the conclusion Fa/Eo 1 as T ^ 0. Thus, we finally conclude that 
both CR and SMS have identical energies Eq at absolute zero in the sense that 
the energy per particle is the same in both states. However, while the entropy of 
CR vanishes, that of the mathematically continued SMS has a negative entropy 
at T = 0. Thus, they are not identical. 

Corollary 4 The free energy Fdis{T), mathematically continued toT — 0, must 
have a maximum at the Kauzmann point T — Tk at which the entropy vanishes. 

From Theorem 1, we know that i^dis(r)=i?o at T = Teq > 0; see point O in 
Fig. 1. From Theorem 3, we also know that Fdis{T)=Eo at T = 0. Because of 
the non-negative heat capacity [d^F^is/dT^ < 0], Fdis(T) is a concave function of 
T. Thus, the mathematically continued -FdisC^) must have a maximum between 
the range (0,Toq) as shown in the inset in Fig. 1. This maximum at T = Tk 
corresponds to the Kauzmann point at which the entropy vanishes. 

Note from Fig. 1, see OO', that the equality Fdis(r) = E'q at T = Toq 
is independent of the way 5'dis approaches zero at i?K- In particular, it does 
not depend on whether there is a singularity in 5dis at £'k, or whether S'dis 
approaches i?K with a finite slope. Thus, the existence of a maximum in Fdis{T) 
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at T = Tk is not a consequence of Theorem 1 alone. We need the additional 
result of Theorem 3. 

The SMS we have defined mathematically between T = and T — Tk 
cannot represent a state of any real system because of the negative entropy, and 
must be replaced by the ideal glass, which is a state that has zero entropy and 
a constant energy Ek and represents the stationary limit of the glassy states 
observed in experiments. Consequently, this stationary limit of the glassy state 
is different not only from the CR, but also from the SMS near absolute zero. 
As SMS is cooled to T = Tk, where it has the energy Ek > Eq, it turns into 
the ideal glass. 

This completes the proof of the generality of our proposed mechanism. In 
CSM, S represents the classical configurational entropy, and in QSM, it rep- 
resents the total entropy Thus, Tk corresponds to vanishing of different 
entropies in the two cases. Any attempt to estimate the classical configurational 
entropy in QSM, where it has no meaning, will require some sort of approxima- 
tion, which we do not consider here. 

The discussion also establishes that the Kauzmann point in the disordered 
phase exists only when there exists another equilibrium state; otherwise, there 
will be no partitioning and, therefore, there will be no Kauzmann point. This 
is most clearly seen in the first example given below. 

6 Exact Model Calculations 

We consider two CSM models in which we obtain positive Kauzmann temper- 
ature. The calculations are carried out exactly. The first example also shows 
how p is used to distinguish different phases, while the second example shows 
that frustration is not necessary for the glass transition. 

6.1 One-dimensional Axis Spin Model 

We now consider a one-dimensional axis spin model, which contains m-component 
spins Si located at site i of the one-dimensional lattice of N sites, with peri- 
odic boundary condition (Sat+i =Si). Each spin can point along or against 
the axes (labeled 1 < A: < m) of an m-dimensional spin space and is of length 
^/^m : S = (0, 0, .., ±-\/to, 0, ..0). The spins interact via a ferromagnetic nearest- 
neighbor interaction energy (—J), with K = J/T > 0. The energy of the inter- 
action is given by 

E = — J Si ■ Si+i. 

i=l...N 

The PF is given by 

Z.iK, m)^[-) E -M-PE) = j T. T , (26) 



20 



where the first sum is over the (2m) ^ spin states of the N spins and T = 
exp(i4rS • S') is the transfer matrix between two neighboring spins. The transfer 
matrix has the eigenvalues [x = exp{K'm)] 

Adis = X + 1/x + 2(m — 1), Aord — X l/x, X — X + 1/x — 2, (27) 

that are 1-fold, m-fold, and (to — l)-fold, respectively 1271. 

We follow de Gennes(28i i29 and provide an alternative and very useful 
interpretation of the above spin model in terms of a polymer system, in which 
each polymer has multiple bonds and loops. The valence at each site in a 
polymer must be even. (The presence of a magnetic field will allow odd valencies, 
which we do not consider here.) The high-temperature expansion of the PF, 
which is given by 



describes such a polymer system, with K > 0, and m denoting the activity 
of a bond and the activity for a loop, respectively, and B and L denoting the 
number of bonds and the number of loops, respectively!^. The empty sites 
represent solvent particles. The number of polymers and the number of bonds 
and loops in each polymer are not fixed and vary according to thermodynamics. 
In addition, there is no interaction between polymers, and between polymers 
and solvent particles, so that the polymer system in H28|l is an athermal solution. 
The temperature T of the spin system does not represent the temperature in 
the polymer problem, as is well known ,,28> i29) . As we will see below, small x 
corresponds to high temperatures where the disordered phase is present, and 
large x corresponds to low temperatures where the ordered and possible SMS 
phases are present. Thus, decreasing T amounts to going towards the region 
where the ordered and metastable disordered phases are present. Let lo denote 
the limiting value as — s- cx) of 



where we have added an uninteresting constant to get rid of the prefactor in (|26|l . 
This is done because the number of microstates appears within the summation 
in the spin model PF in H26(l . Thus, the inclusion of the prefactor will make 
the microstate entropy negative. The prefactor is, however, required for the 
polymer mapping. 

The importance of the polymer mapping is that we can take to > to be 
a real number, even though non-integer to makes no sense for a physical spin. 
Thus, for non-integer values of to, only the polymer system represents a physical 
system. For m — 1, the axis model reduces to the Ising model, while for to — > 0, 
it reduces to the a model of linear chains with no loops |28[ 12^ . The eigenvalue 
Adis is dominant at high temperatures for all to > and describes the disordered 
phase. Its eigenvector is 




(28) 



UJ = [l/N) lnZjv(A',m) + ln(2TO), 



(29) 
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where {2k\ ( or {2k + 1|) denotes the single-spin state in which the spin points 
along the positive (or negative) k-th spin-axis. It has the correct symmetry 
to give zero magnetization (p — 0). For m > 1, Adis remains the dominant 
eigenvalue at all temperatures T > 0. For < m < 1, the situation changes and 
Aord becomes dominant at low temperatures T < Tc, or [x > Xc = 1/(1 — m)] 
where Tc is determined by the critical value Xc = exp(Jm/Tc); there is a phase 
transition at Tc- The corresponding eigenvectors are given by the combinations 

= [(2A:|-(2A:+1|]/V2, = 0, 2, .., m - 1, 

which are orthogonal to (^dis | , as can be easily checked. These eigenvectors have 
the symmetry to ensure p ^ 0. The remaining eigenvalue A is (m — l)-fold 
degenerate with eigenvectors 

(J^C^+i) I = [(2fc| + (2fc + 1| - {{2k + 2\- {2k + 3|)]/V4, k^0,2,..,m- 2. 

For m > 0, this eigenvalue is never dominant. For m — *■ 0, it becomes degenerate 
with Adis. Since the degeneracy plays no role in the thermodynamic limit, there 
is no need to consider this eigenvalue separately for m > 0. 

We now consider the limit iV — > oo. The adimensional free energy per site, 
which represents the osmotic pressure |S1 EO], of the high-temperature equilib- 
rium phase is WdisC?") = In(Adis)- It can be continued all the way down to T = 0, 
even though the equilibrium osmotic pressure has a singularity at Xc- Similarly, 
'A'ord(T) = In(Aord) related to the low-temperature equilibrium phase can be con- 
tinued all the way up to T ^ oo. To calculate the entropy density, we proceed 
as follows. The bond and loop densities are given by 

0B = duj/d In K, (j)L = duj/d In m, (30) 

which are needed to calculate the entropy per site of the polymer system 

s^^-* = cj — 03 In — 0L In m; 

the superscript is to indicate that it is the polymer system entropy, and is 
different from the spin system entropy s^^^ — dTuj/dT. If we define uj without 
the last term in H29|) . then 0l and s^^-* must be replaced by — 1) and 
(s^^) — In 2), respectively. This will not affect any of the conclusions below. 

In the following, we will be only interested in the polymer entropy. The 
proper stability requirements for the polymer system are 

{d(t>B/d\nK) > 0, (a0L/91nm) > 0, (31) 

as can easily be seen from H26(l . and must be satisfied even for SMS. They replace 
the positivity of the heat capacity of the spin system, which no longer represents 
a physical spin system for < m < 1. It is easy to see from the definition of 
°dis t^^^t (^^dis /^^)ni need not be positive, even if the conditions in l|31|) are 
satisfied. 
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It should be noted that the Theorem 3 was for a canonical PF, whereas 
the polymer PF in H28() is an athermal grand canonical PF. The proof of the 
theorem can be easily extended to this or any other ensemble with similar results. 
Here, we will instead give a direct demonstration of the theorem. For this, we 
compute Lu as K oo {T 0) for the two eigenvalues Adis and Aord- From 
H3U|) . it is easy to see that (j>B ^ mK for both states as T — > 0. Thus, using 
w = s'^' + (/)B In K + In m, we have 

Wdis(r)/tJord(T) ^ 1 as T^O. (32) 

This is in accordance with H25I) . This means that if the eigenvalue Adis is taken 
to represent the metastable phase above Xc, its osmotic pressure must become 
equal to that of the equilibrium phase (described by the eigenvalue Aord) at 
absolute zero, in conformity with Theorem 3. We take Wdis(T') to represent the 

(S) 

SMS osmotic pressure below T^. One can also check that Ts^;,' ^ 0, as T ^ 0. 

We will only discuss the disordered polymer phase below for < m < 1. It 
is easily checked that the above stability conditions in (|31|l are always satisfied 
for Adis! see, for example, the behavior of in Fig. 2, where we have taken 
m = 0.7, and J — \. Since the high-temperature disordered phase represents a 
physical system, it cannot give rise to a negative entropy s^jj^ above Tc\ how- 
ever, its metastable extension violates it as shown in Fig. 2, where its entropy 
(p) 

Sjjg becomes negative below Tk = 0.266, which is lower than the transition 
temperature Tc. 

We now make an important observation. As m decreases (below 1), both 
Tk and Tc (Tk < Tc) move down towards zero simultaneously. As m ^ 0, the 
equilibrium ordered phase corresponding to Aord disappear completely, and the 
disordered phase corresponding to Adis becomes the equilibrium phase. There 
is no transition to any other state. Thus, there is no metastability anymore. 
Consequently, there is no Kauzmann point since there is no other ordered state 
any more, as argued above. Thus, our exact calculation confirms our earlier 
conclusion that the existence of an ordered state is crucial for the existence of the 
entropy crisis. The existence of an ordered state sets the zero of the temperature 
scale by its minimum energy Eq. This scale then sets the temperature Tk of 
the lowest SMS energy > Eq to be positive. Thus, one must consider the 
ordered and the metastable states together. 

We also observe that there is no singularity in Adis or LJtiXs{T) at Tc, even 
though there is a phase transition there. Similarly, there is no singularity in Aord 
or Word(T') at Tc. Thus, the thermodynamic singularity in the equilibrium free 
energy does not necessarily create a singularity in lj^{s{T) or Word (7") at Tc, as 
was discussed earlier. The existence of a singularity or spinodal at some other 
temperature is a different matter. 

It should be noted that the eigenvalues Adis and Aord are independent of 
the size of the lattice. Therefore, they can be used to describe not only the 
disordered and ordered phases, but also the SMS, which is the continuation of 
the disordered phase, even for a finite N . Thermodynamic limit is not necessary. 
For finite N, Za{K,m) < Z(K,m) and the inequality becomes an equality in 
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the sense of (@J) only as A?^ oo for the proper choice of a depending on the 
temperature. 

6.2 Binary Mixture Model 

We now consider a simple lattice model of an incompressible binary mixture of 
two kinds of particles A and B, to be represented by an Ising spin S. The two 
spin states (+1 or up) and (—1 or down) represent the particles of two species 
A and B, respectively. As we are not interested in their phase separation, but in 
the possibility of a glass transition, we assume that their mutual interaction is 
attractive. In addition, we are interested in a first-order transition for conven- 
tional supercooling. We will, therefore, use an anti-ferromagnetic Ising model 
in zero magnetic field with both two-spin (J > 0), and three-spin interactions 
(J' 7^ 0); the latter ensures that the melting transition is first order. In order 
to solve the model exactly, we consider a Husimi cactus made of squares, on 
which the model can be solved exactly We consider the simplest cactus in 
which only two squares meet at a site; they cannot share a lattice bond. The 
squares are connected so that there are no closed loops except those formed by 
the squares. The cactus can be thought as an approximation of a square lattice, 
so that the exact Husimi cactus solution can be thought of as an approximate 
solution of the square lattice model. There is a sublattice structure at low tem- 
peratures caused by the anti-ferromagnetic interaction: particles of one species 
are found on one of the two sublattices. We identify this ordered structure as a 
crystal. The interaction energy is 

E^J^ SS' + J'Y, ■ (33) 

The first sum is over nearest-neighbor spin pairs and the second over neighboring 
spin triplets within each square. In the absence of the three-spin coupling, 
the two-spin coupling gives rise to an antiferromagnetic (AF) ordering at low 
temperatures. For J' > (—J), the AF ordering remains the preferred ordering, 
while for J' < {—J), the ferromagnetic ordering is prefereed. Threrefore, we 
only consider J' > {—J) in the following. We set J=l to set the temperature 
scale. 

The model is solved recursively, as has been described elsewhere |2l-{j. We 
label sites on the cactus by an index m, which increases sequentially outwards 
from m = at the origin. We introduce partial PF's Z„i{1) and Zm{i), de- 
pending on the states of the spin at the m-th cactus level. It represents the 
contribution of the part of the cactus above that level to the total PF. We then 
introduce the ratio 

(T)/[Z,„(T)+^m(i)], (34) 

which satisfies the recursion relation 

Xyn = f{Xm+l,Xm+2,v)/[S{Xyn+l,Xyn+2,v) /(y^+i , ?/m+2 , l/w )] , (35) 
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where 



f{x, x' , v) = x^x' /u^v^ + 2xx'yv^ + x^yv^ + u^x'y"^ + 2xyy' + y'^y' /v^ , (36) 
with 

It = , V = e^"^ ,y = 1 ~ x,y' = 1 — x' . 

There are two kinds of fix-point solutions of the recursion relation that de- 
scribe the bulk behavior |23|. In the 1-cycle solution, the fix-point solution be- 
comes independent of the level index m as we move towards the origin m — 0, 
and is represented by x* . For the current problem, it is given by x* — 1/2, as 
can be checked explicitly by the above recursion relation in 135|) . It is obvi- 
ous that it exists at all temperatures. There is no singularity in this fix-point 
solution. This solution corresponds to the disordered paramagnetic phase at 
high temperatures and the SMS below the melting transition. The other fix- 
point solution is a 2-cycle solution, which has been found and discussed earlier 
in the semi-flexible polymer problem[71 the dimer model|3J, and star and 
dendrimer solutions [32] ■ The fix-point solution alternates between two values 
xl, and X2 on two successive levels. At T = 0, this solution is given either 
by xl = 1, and = 0, or by = 0, and = 1. The system picks one of 
these as the solution. At and near T = 0, this solution corresponds to the low 
temperature AF-ordered phase, which represents the CR and its excitation at 
equal occupation, and can be obtained numerically. The 1-cycle free energy is 
calculated by the general method proposed in and the 2-cycle free energy 
is calculated by the method given in [711^. 

For J' — 0.01, we find that Tm = 2.753, where there is a discontinuity in the 
entropy per site of 0.0168. The SMS below Tm represents SCL, whose entropy 
density, see Fig. 2, vanishes at Tk = 1.132, and whose specific heat (not shown) 
remains positive everywhere with a maximum at T = 1.26. At absolute zero, 
the entropy per site Sdis — —0.3466, while the CR entropy is zero, as expected. 
The CR and SCL free energies per site become identical {— 2 J) at absolute zero 
in accordance with the Theorem 3. Thus, the free energy diagram we obtain in 
this case is similar to that in the inset in Fig. 1. 

7 Discussion &; Conclusions 

7.1 Thermodynamic Criterion for Ideal Glass Transition 

The work was motivated by a desire to identify a general thermodynamic mech- 
anism for the glass transition. For this, we identify a general thermodynamic 
condition. This general principle has also been verified in some recent work 
on lattice models that have been carried out|31[ 15^. and has its foundation 
in the original idea of the entropy crisis S'ex < noted by Kauzmann, and 
is as follows. The entropy due to a set of coupled degrees of freedom, when 
properly defined either using quantum mechanics or by discretization as in a 
lattice model, can never be negative, since the number of configurations for a 
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given set of macroscopic quantities must be an integer > 1. Consequently, this 
condition of non-negative entropy S is nothing but the condition of reahty. Its 
violation gives rise to the concept of absolute entropy crisis. It is in this sense 
we have used the entropy crisis in this work. Such a notion is more stringent 
than the Kauzmann criterion that Sex be not negative, for which there is no 
thermodynamic justification, as we have argued in the Introduction. For the 
glass transition, we are interested in the set of coupled degrees of freedom that 
contains the configurational (i.e., positional) degrees of freedom. The entropy 
of such a set, regardless of whether quantum or classical mechanics is used in 
its calculation, is what we call the configurational entropy. 

It should be noted that there are various other definitions of the configura- 
tional entropy in the literature. Many workers take Sex to denote the configura- 
tional entropy. However, as said above, there is no thermodynamic basis for Sex 
to be non-negative. So, it cannot play any role in a thermodynamic theory of 
glass transition. The landscape picture identifies the entropy 5*18 of the inherent 
structure with the configurational entropy. Its vanishing is used to identify the 
glass transition. It is easy to see|33] that for the classical configurational PF, 
the two entropies are related: 

= 5is(r) + 5basi„(T), 

where S'basin(7") = 9(rinZbasin)/9T is the entropy arising from the average 
basin PF .^basin(2^) in the landscape picture. Thus, our criterion S{T) = of 
the entropy crisis is also stringent than Sis{T) = in the landscape picture. Our 
criterion will also require Shasin{T) = 0, which can only occur at a temperature 
lower than the temperature at which Sis{T) — 0. It should be noted that there 
is no kinetic energy contribution in 5'basin(T') as the landscape picture deals only 
with the classical configurational PF; the translational degrees of freedom are 
decoupled in classical mechanics, as discussed above. 

7.2 Continuum vs. Discrete Models 

The best known example of classical models giving rise to negative entropy at 
low temperatures is the ideal gas. Similarly, classical real gases also give rise 
to negative entropy at low temperatures. This problem can be easily traced to 
the fact that we are treating the real and momentum spaces as continuum jUll. 
Another well-known example is the Tonks gas of rods in one-dimension (no 
kinetic energy), which also gives negative entropy at high coverage [TTll31| . Here, 
the one-dimensional space is treated as continuum. On the other hand, it is 
well known that a classical lattice model will never give rise to negative entropy. 
Similarly, the random energy model 12 , which treats energy as continuum, gives 
negative entropy at low temperatures. 

It is clear that the problem of negative entropy is not due to the classical 
nature of the systems, but due to the continuum nature of the model|TS|. To 
ensure non-negative entropy, we must discretizc the model, as we have discussed 
here. Once this has been done, the entropy crisis becomes a genuine crisis 
imposed by the reality condition as we have proposed earlier 
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7.3 General Thermodynamic Mechanism for Ideal Glass 
Transition 



The general thermodynamic mechanism of the ideal glass transition occurs in 
any system that has an ordered state, distinct from the high-temperature dis- 
ordered state. The continuation of the free energy of the disordered state below 
the melting transition at Tm gives the free energy of the stationary metastable 
state. We have shown that this continuation always gives rise to a stable free 
energy. For example, it gives rise to a non-negative heat capacity for T > 0. 
This SMS free energy fdis{T) becomes equal to Eq at T^q > 0, as the temper- 
ature is reduced. The energy continues to drop until finally, fdis{T) = Eq at 
T — 0. Since /dis(r) — Eq at T — and at T = Tcq, it must have a maximum 
at some intermediate temperature < Tk < 2eq- The energy Ek > Eq at Tk 
because of positive heat capacity. The entropy is zero at Tk and negative below 
Tk. Thus, the SMS over the range (0,Tk) is unphysical and must be replaced 
by the ideal glass. 

The above mechanism has been shown to be generic by the rigorous analysis, 
which is valid for classical and quantum systems. To the best of our knowledge, 
this result is the first of its kind and shows that the entropy crisis is genuine in 
those systems in which there is a more stable ordered state than the disordered 
SMS. 

We have also considered models in which there are two possible transi- 
tions; one of them is from a disordered phase to an intermediate phase, and the 
second one at a lower temperature from the intermediate phase to an ordered 
phase. In this case, two possible SMS's as the continuation of the disordered 
and intermediate phases emerge, with, each extension giving rise to its own 
entropy crisis; in addition, both have the same free energy as T = as the 
ordered phase. Thus, the mechanism is generic. 

7.4 SMS & Exact Calculations 

The transition between SMS and the ideal glass is not brought about by any 
thermodynamic singularity at Tk; rather, it is imposed by the reality require- 
ment. To the best of our knowledge, the ideal glass state does not explicitly 
emerge as a phase in any calculation that has been carried out so far, includ- 
ing the two that are presented here. In this sense, this transition is a very 
special kind of transition, which does not seem to belong to the class of phase 
transitions in which various phases emerge in the calculation. 

The two examples that we have presented here show the existence of SMS. 
Thus, they demonstrate that our hypothesis of SMS existence is not vacuous. 
Both examples also show genuine entropy crisis in SMS. Thus, they provide 
support for the violation of the reality condition in SMS at a positive temper- 
ature in exact calculations. We need to invoke an ideal glass transition at this 
temperature in each model. 

The one-dimensional exact calculation is not a mean-field type calculation, 
and is presented not only to overcome the folklore that SMS's exist only at the 
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mean-field level, but also to explicitly confirm the theorems. It is a model of a 
polymer system and confirms all the theorems. In particular, it shows that as 
m — > 0, the ideal glass transition disappears because there is no ordered state 
anymore. Hence, our results are not of mean-field nature. 

The second example is also presented not only to overcome the folklore that 
frustration is crucial in the transition, but also to establish that entropy crisis 
is possible in systems containing small molecules, and not just polymers. There 
is no frustration in the Ising model because the cactus consists of squares. This 
example also confirms the theorems and the corollary. 

Both examples show that the ideal glass transition occurs at a positive tem- 
perature and the energy Ek at that temperature is higher than Eq, the ideal 
CR energy at T = 0. This means that the ideal glass has a higher energy than 
the crystal at absolute zero, in conformity with the experiments. 

7.5 Absence of Entropy Crisis in Equilibrium State 

We now come to a very important consequence of our corollary. It is the fol- 
lowing. The presence of an entropy crisis at a positive temperature implies that 
there must exist an equilibrium state for which no entropy crisis can exist. The 
equilibrium states in any system or model calculation, if carried out exactly, 
will never give rise to any entropy crisis at a positive temperature. This is 
because the lowest energy Eq determines the lowest allowed temperature T = 
in the system, even for a finite system. Since there cannot be any singularity 
in any finite-system equilibrium free energy, the latter should continue all the 
way down to T = and the equilibrium entropy must remain non-negative at 
all temperatures; the latter can only vanish at T = 0. Thus, there would be 
no entropy crisis at a positive temperature in the equilibrium state. This will 
remain true as the thermodynamic limit is taken. Thus, no entropy crisis can 
occur at a finite temperature in the equilibrium state in that the entropy be- 
comes negative below that temperature. However, it is possible that as N oo, 
the equilibrium free energy becomes horizontal, so that the entropy vanishes, 
over a non-zero temperature range (0,Tc). Since the equilibrium free energy 
exists with non-negative entropy for all T > for finite N, the free energy must 
show a singularity at Tq- While the system is frozen over the range (0, Tq), its 
appearance is accompanied by a phase transition. (This should be contrasted 
with the existence of a Kauzmann point, below which the entropy becomes neg- 
ative, but its appearance is not accompanied by any singularity in the SMS free 
energy. Replacing the unphysical SMS free energy below the Kauzmann point 
by a frozen state is done by hand; it docs not emerge as part of the calcula- 
tion.) Thus, we conclude that equilibrium state in any system will never show 
an entropy crisis at a positive temperature. The zero of the temperature scale 
is determined by the lowest possible energy Eq. If any exact calculation for the 
free energy or the entropy predicts an entropy crisis at a positive temperature, 
this will necessarily imply that there must exist another state, the equilibrium 
state, which will not show an entropy crisis. 

This observation has been crucial in a recent investigation of a dimer model 
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|5T] in which the disordered phase underwent a first-order transition to an equi- 
librium ordered phase. The ordered phase then gave rise to an entropy crisis 
at a lower temperature, which forced us to look for another equilibrium state, 
which was eventually discovered above the temperature where the entropy crisis 
was found, so that the crisis occurred in a metastable state (this time emerg- 
ing form an intermediate ordered state) as we have suggested and followed the 
mechanism proposed here. 

7.6 Landscape Picture 

Finally, we wish to make connection of the ideal glass energy i?K with the 
inherent structure in the landscape picture. The ideal glass at T = in the 
canonical ensemble must be at a local minimum in the landscape. Hence, its 
energy must be the energy of the particular inherent structure. Since the ideal 
glass emerges at Tk, where S{T) — 0, we are forced to conclude that this 
inherent structure also represents SMS at Tk. This most certainly implies that 
SMS must be confined in the basin associated with the SMS inherent structure. 
This confinement must occur at a higher temperature (than Tk) where Sis(T) = 
0. 

In conclusion, we have justified the mechanism that gives rise to an entropy 
crisis in metastable states in systems in which there exists is a more stable phase. 
The generality of the mechanism is refiected in the generality of the validation of 
the mechanism, which is common in classical and quantum mechanical systems. 

We would like to thank Andrea Corsi and Fedor Semerianov for various 
useful discussions, and help with the figures (Andrea Corsi). 
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Figure Captions 

1. Schematic form of the generic entropy functions for various possible states. 

2. The bond and the entropy densities. Both models show an entropy crisis 
at a positive temperature. 
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